Geometrical scaling in hadronic collisions* 
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We show that the pt spectra measured in pp collisions at the LHC 
exhibit geometrical scaling introduced earlier in the context of deep inelastic 
scattering. We also argue that the onset of geometrical scaling can be seen 
in nucleus-nucleus collisions at lower RHIC energies. 



1. Introduction 

With the start of the LHC we have been confronted with a wealth of data 
on multiparticle production at high energies both in pp PQ-[3] and in heavy 
ion collisions [1] . One of the remarkable results is that total multiplicity of 
charged particles produced in central rapidity in pp collisions is rising like 
a power of s 

~ s x with A ~ 0.23. (1) 

The power law behavior ([I]) is expected in saturation models [5]-|10|. In 
this paper we show that power like behavior of total multiplicity follows 
naturally if pi spectra of charged particles exhibit geometrical scaling 
E2]. In Sect. [| we remind basic properties of geometrical scaling which 
was introduced in the context of small x deep inelastic scattering (DIS). 
In Sect. [3] we introduce geometrical scaling in pp collisions. In Sect. [4] we 
briefly discuss a possibility of geometrical scaling in heavy ion collisions. 
Finally we summarize and give conclusions in Sect. [5] 

* Based on a talk presented at Cracow Epiphany Conference " On the First Year of the 

LHC", Krakow, January 10-12, 2011. 
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2. Geometrical scaling in DIS 

In a successful description of small x DIS proposed in seminal papers by 
Golec-Biernat and Wiisthoff (GBW model) [13], a cross-section for virtual 
photon-proton scattering in DIS reads: 

a rp = J dr 2 \^{r,Q 2 )\ 2 a dP (r 2 Ql(x)). (2) 

Here ip is the wave function describing dissociation of a virtual photon into 
a qq dipole and a^p is a dipole-proton cross-section. The main assumption 
of the GBW model is that adp which in principle is a function of two 
independent variables: dipole size r and dipole-proton energy W (or Bjorken 
x), depends in practice only on a certain combination of these two variables, 
namely on the product r 2 Q 2 (x) where 

Q 2 s (x) = Q 2 (^) X (3) 
is called a saturation scale. Bjorken x is defined as 

X = Q^+W 2 ' (4) 

Here Qo ~ 1 GeV and xq ~ 10 -3 are free parameters whose precise values 
can be extracted by fitting ^ to the HERA data. Power A is known to be 
of the order A ~ 0.2 -=- 0.3. 

For transverse photons (neglecting quark masses): 

\ih(r, Q 2 )\ 2 = T dz [z 2 + (1 - z) 2 ] Q 2 K 2 (Qr) (5) 
Jo 

where 

Q 2 = z{\ - z)Q 2 (6) 
and K\ is a modified Bessel function. From Eq. ([5j it follows that 

\Mr,Q 2 )\ 2 = Q 2 \MrQ)\ 2 (7) 

where we have explicitly factored out Q 2 . Defining new variable u = Q 2 r 2 , 
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new function 



and scaling variable r 
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we arrive at: 

o- 7 *p = J -^<f)(u)a d p(uT). (9) 

It follows that (7-y*p IS £1 function of scaling variable r, rather than a function 
of two variables Q 2 and x. This phenomenon is known as geometrical scaling 
(GS)[H]. GS has been observed in DIS data for x < 0.01. 
In the Golec-Biernat-Wiisthoff model 

<?dp{r 2 Ql(x)) = a (1 - exp(-r 2 Q s 2 (x))) (10) 

where ctq is dimensional constant; <7o — 23 mb. 




Fig. 1. Dependence of A c ff on Q 2 from HERA (HERA data points [H] after 
Ref.QT]). 

In practice Q 2 may also have some residual dependence on Q 2 if DGLAP 
evolution in ^ is taken into account [15]. Indeed, it can be shown that in 
the GBW model - up to logarithmic corrections - DIS structure function is 
proportional to Q 2 : 

a rp (x, Q 2 ) ~ ^o^^ and F 2 [x, Q 2 ) ~ a Q 2 (x). (11) 



At first sight Eq. (11) may look contradictory, since left hand sides depend 
on Q 2 and the right hand sides do not. In practice, exact calculation of the 
integral in ^ renders some mild Q 2 dependence of the right hand sides. 
Moreover DGLAP evolution introduces Q 2 dependence of Q 2 (x). Therefore 
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the effective saturation scale can be conveniently parameterized as 

qU*.q 2 ) = qU^Y mq ' } . (12) 



X 

}2 



Exponent A e ff(Q ) has been extracted from the HERA data [16J. This is 
shown in Fig. [T] (after Ref.|17j) together with an eyeballing fit 



/o 2 \ - 35 

A eff (Q) = 0.13 + 0.1 f^J . (13) 

3. Geometrical scaling in pp collisions 

In pp collisions we do not have a viable model of low and medium p^ 
particle production at high energies. Nevertheless one often uses a fcj factor- 
ized form of a cross-section describing production of a pt gluon at rapidity 

d(J 37T 1 f 

E d?p = Yp^ J dk T a s(k T )^ 1 (x 1 ,kl)ip 2 (x2, (k ~p)t) (14) 

where 

xi,2 = ^e±y (15) 

are Bjorken x's of colliding partons. Here 92's are unintegrated gluon densi- 
ties. Introducing "regular" gluon distribution 

Q 2 

xG(x,Q 2 ) = J ' dk^(x,k\). (16) 



one obtains for p\ > Q 2 



E ~¥p = i " 5 ^ {fi{xi,p\)x 2 G(x2,p%) + Lp 2 {x 2 ,p\)xiG{xi,p\)} . 
There have been recently more involved model calculations of particle mul- 



tiplicity based on Eq. ([14|) [T9]-[2TJ. 

Kharzeev an Levin proposed a simple Ansatz for unintegrated gluon 
distribution [22J: 



3a f 1 for p T < Ql 
vras(Q s ) 1 Q 2 /p 2 T for Q 2 <p 2. 
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Hence, up to the logarithmic corrections due to the running coupling con- 
stant, we arrive at geometrical scaling for the multiplicity distribution 

^ 1 E%-=*F{t) (19) 



drjd 2 pT cr ine i d 3 p Ql 
where Qo ~ 1 GeV and o"i ne i is the inelastic cross-section. Although we 



have used a very simple Ansatz (18) for unintegrated gluon distribution tp, 
it satisfies the generic property that xG(x, ~ Q 2 which is enough for 
GS to hold. 




Fig. 2. Charged particle multiplicity at mid rapidity < 2.4 as measured by CMS 
0, plotted as functions of p\ (left) and scaling variable r (right) for A = 0.27. 



(20) 



F(t) is a universal function of the scaling variable 

T= v\ 
Q 

where in view of §5§ and ([15]) 

where W = yfs x 10 -3 . Here factor 10 -3 corresponds to the (arbitrary at 
this moment) choice of xq. 

The power like growth of the multiplicity can be easily understood as a 
consequence of geometrical scaling. Indeed 



dy J Ql 



2 F(r). (22) 
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Simple change of variables gives [TT] 

dp\ _ 2 fW 



2A 

2+A A 

r 2+A(ir. (23) 



The integral over dr is convergent and universal, i.e. it does not depend on 



energy. It follows from Eq. ( 23 ) that the effective power of the multiplicity 
growth is 

A=^<A (24) 

rather than A. For A = 0.27 we have that A = 0.238. 

In Refs.[llJ it was shown that CMS charged particle pp spectra [2] at 
mid rapidity \rj\ < 2.4 plotted as functions of scaling variable r fall on one 
universal curve (19). This is depicted in Fig. [2] where we plot pr spectra 



for three LHC energies as functions of p\ (left panel) and as functions of 
scaling variable r for A = 0.27 (right panel). 

In order to examine the quality of geometrical scaling in pp collisions 
we plot in Fig. [3] ratios of spectra measured at 7 TeV to spectra at 0.9 and 
2.36 TeV in function of (left panel) and \fr (right panel). We see that 
original ratios plotted in terms of pt range from 1.5 to 7, whereas plotted 
in terms of yfr they are well concentrated around unity. This is further 
illustrated in the left panel of Fig. [4] which presents the enlarged view of the 
right panel of Fig. [3j With this accuracy we see a small systematic increase 
of the ratios (apart from the first 4, 5 points) which suggests some weak 
dependence of exponent A on pp- The value of A = 0.27 has been obtained 
by minimizing deviations of ratios -R7/0.9 an d -R7/2.36 from 1 for central pr 
points [i.e. rejecting first 5 and last 4 points). 

Residual dependence of exponent A on pi is in agreement with small x 
dependence of the DIS structure function as measured in HERA [16] . In 
Ref. [12] we have argued that this dependence can be well approximated 
by use of the effective exponent A e g- of Eq. ( p~3| ) with argument Q = 2pp- 
This is demonstrated in the right panel of Fig. [4] where we used A e ff(2pr) 
to calculate the ratios -R7/0.9 and -^7/2.36- An impressive improvement of 
geometrical scaling (i.e. of the equalities -R7/0.9 — 1 and -R7/2.36 — 1) ca n 
be indeed seen. 



4. Onset of geometrical scaling in heavy ion collisions 

Heavy ions provide much readier information on the characteristics of 
particle production at high energies. Indeed, one can study not only energy 
dependence but also atomic number A— dependence, rapidity dependence 
(at RHIC much larger rapidity range has been covered than so far at the 
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Fig. 3. Ratios of CMS p T spectra [2] at 0.7 TeV to 0.9 (blue circles) and 2.36 TeV 
(red triangles) plotted as functions of pt (left) and scaling variable y/r (right) for 
A = 0.27. 





A = 3. Bls (2p T ) 



Fig. 4. Enlarged plot of the right panel of Fig. [3] for A = 0.27 (left) and for A 

Aeff(2p T )- 



LHC) and finally centrality dependence. The production of quark-gluon 
plasma and its ability to "remember" the initial conditions of the saturated 
gluonic matter are here of primary interest. Unfortunately RHIC energies 
are presumably too low for geometrical scaling to work. Nevertheless we 
show below, that approximate GS can be seen in the RHIC data. To this 
end we choose the PHOBOS px distributions measured in gold-gold and 
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copper-copper collisions at 62.4 and 200 GeV per nucleon 

Here a new scaling law is particularly interesting. Namely the satura- 
tion scale in nucleus-nucleus collisions scales with A as [25] (for review see 
Ref. [26]): 

Q\ s = A l '*Ql (25) 
which implies that the relevant scaling variable reads: 

A ~ AV3Q2 - A 1/3Q2 [ W ) ■ W 

In Fig. [5] we plot multiplicity distribution for central Au-Au and Cu-Cu 
collisions in function of p\ and ta- In this case a slightly higher value of 
the exponent A is used, namely A = 0.3. We see again that the rescaled 
spectra seem to fall on one curve, although the alinement of Au and Cu 
spectra is not perfect for small and medium values of ta- Nevertheless a 
tendency towards geometrical scaling is clearly seen. Similar conclusions 
can be drawn for more peripheral collisions. A detailed study of the onset 
of geometrical scaling in heavy ion collisions will be presented elsewhere. 



PHOBOS - 6 % 
■ Au 200 GeV 

□ Au 62.4 GeV 
A Cu 200 GeV 

A Cu 62.4 GeV 
A =0 




PHOBOS - 6 % 
■ Au 200 GeV 

□ Au 62.4 GeV 
A Cu 200 GeV 
A Cu 62.4 GeV 
A = 0.30 



PrlGeV/c] 



Fig. 5. Multiplicity distribution in heavy ion collisions for Au-Au and Cu-Cu at 
two RHIC energies 200 and 62.4 GeV EH CU] plotted in terms of p\ (left panel) 
and scaling variable (right panel). 



5. Conclusions 



In this paper we have demonstrated that geometrical scaling originally 
postulated in deep inelastic scattering [T3] is also exhibited by the px spectra 
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in hadronic collisions |11] [T2] . To this end recent CMS data [2] have been 
analyzed and shown to scale with scaling variable r defined in Eqs.( 
simplified model of Gribov, Levin and Ryskin [18J has been used to motivate 
the appearance of GS in hadronic collisions. This model can be a'priori used 
to study the shape of the universal scaling function F(r) which deserves a 
separate study. 

A notable difference between DIS and hadronic collisions is that in DIS 
we deal with totally inclusive cross-section, whereas in pp both hadroniza- 
tion and final state interactions play essential role. Nevertheless the imprint 
of the saturation scale Q s is visible in the spectra, which means that the 
information on the initial fireball survives until final hadrons are formed. 

It has been shown that the quality of geometrical scaling is improved if 
the exponent A becomes pi-dependent [12] in accordance with Q-dependence 
of \ e s(Q = 2pr) obtained from DIS. This is a remarkable feature that sup- 
ports the picture in which medium particles are produced from saturated 
gluonic matter irrespectively of the scattering states. 

If so, geometrical scaling should be also present in heavy ion collisions. 
The detailed studies will be certainly carried out at the LHC. Here we 
have analyzed PHOBOS data |23} [23] for two RHIC energies and for two 
different nuclei: gold and copper, and the onset of geometrical scaling has 
been clearly seen. Interestingly, we have found that the exponent A that 
governs geometrical scaling in heavy ion collisions is higher than the one 
in pp. This is in striking agreement with the fact that multiplicity growth 
with energy observed by ALICE [3] is faster in heavy ions than in pp. 
Question arises to what extent the hydrodynamical evolution of the quark- 
gluon plasma is going to wash out geometrical scaling that is present in 
the initial state. Further studies should also concentrate on centrality and 
rapidity dependence of GS. 
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